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Abstract — Integer-forcing (IF) linear receiver has been recently 
introduced for multiple-input multiple-output (MIMO) fading 
channels. The receiver has to compute an integer linear com- 
bination of the symbols as a part of the decoding process. 
In particular, the integer coefficients have to be chosen based 
on the channel realizations, and the choice of such coefficients 
is known to determine the receiver performance. The original 
known solution of finding these integers was based on exhaus- 
tive search. A practical algorithm based on Hermite-Korkine- 
Zolotareff (HKZ) and Minkowski lattice reduction algorithms 
was also proposed recently. In this paper, we propose a low- 
complexity method based on complex LLL algorithm to obtain 
the integer coefficients for the IF receiver. For the 2x2 MIMO 
channel, we study the effectiveness of the proposed method in 
terms of the ergodic rate. We also compare the bit error rate 
(BER) of our approach with that of other linear receivers, and 
show that the suggested algorithm outperforms the minimum 
mean square estimator (MMSE) and zero-forcing (ZF) linear 
receivers, but trades-off error performance for complexity in 
comparison with the IF receiver based on exhaustive search or 
on HKZ and Minkowski lattice reduction algorithms. 

Index Terms — CLLL algorithm, MIMO, linear receivers. 

I. Introduction 

In multipath fading channels, using multiple antennas at the 
transceivers is known to provide large capacity gains. Such a 
capacity gain comes at the cost of high decoding complexity 
at the receiver. It is known that a high-complexity joint ML 
decoder can be employed at the receiver to reliably recover 
the information. On the other hand, the linear receivers such 
as the ZF and the MMSE receiver J6) reduce the decoding 
complexity trading-off error performance. 

The integer forcing (IF) linear receiver has been recently 
proposed fl4l . This new architecture obtains high rates in 
MIMO fading channels. In this approach, the transmitter em- 
ploys a layered structure with identical lattice codes for each 
layer. Then each receive antenna is allowed to find an integer 
linear combination of transmitted codewords. The decoded 
point will be another lattice point because any integer linear 
combination of lattice points is another lattice point. This idea 
has been brought to MIMO channels from the compute-and- 
forward protocol for physical layer network coding E], El- 

In the MIMO IF architecture, a filtering matrix B and a 
non-singular integer matrix A are needed such that BH sa A 
with minimum quantization error at high signal-to-noise ratio 
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(SNR) values. The exhaustive search solution to the problem of 
finding A,B is addressed in JT4). It is prohibitively complex 
already for 2x2 real MIMO and becomes untractable for 
2x2 complex MIMO and beyond. A smart practical method 
of finding A based on HKZ and Minkowski lattice reduction 
algorithms has been proposed recently j8]. This provides full 
receive diversity with much lower complexity in comparison 
to exhaustive search. The major differences between integer- 
forcing linear receivers and lattice reduction aided MIMO 
detectors iflOl . ifTTI . ifTSl are also presented in fl8]. 

In this paper, we propose a low-complexity method for 
choosing the above matrices. In J8), [14), a 2n,-layered scheme 
is considered with real lattice codebook for each layer. Unlike 
the model there, we work on complex entries and we lift that 
set-up to complex case. The proposed method is a combination 
of three low-complexity methods which are based on complex 
lattice reduction (CLLL) J2) technique for a lattice, and 
singular value decomposition (SVD) of matrices. For the 2x2 
MIMO channel, we compare the performance (in terms of 
ergodic rate and uncoded probability of error) of the proposed 
low-complexity solution with the known linear receivers and 
show that the proposed solution (;) provides a lower bound 
on the ergodic rate of the IF receiver, (ii) outperforms the ZF 
and MMSE receivers in probability of error, and (Hi) trades off 
error performance for computational complexity in compari son 
with exhaustive search and other lattice reduction methods 
including HKZ and Minkowski algorithms. 

The rest of the paper is organized as follows. In Section [U] 
we give a brief background on lattices. We present the problem 
statement along with the signal model in Section [In] In Sec- 
tion Hy] we study the solution to the IF receiver via two CLLL 
algorithms. A complexity comparison for different known 
approaches is also given in this section. In Section [V] we show 
some simulation results on the performance of IF receiver in 
ergodic MIMO setting. Finally, we present concluding remarks 
in Section [VI] 

Notations. Row vectors are presented by boldface letters, 
and matrices are denoted by capital boldface letters. Let v 
be a vector, v T denotes transposition, and v h denotes the 
Hermitian transposition. We show the nx n identity and zero 
matrix as I„ and 0„ respectively. For a matrix X, the element 
in the k-th row and m-th column of X will be denoted by 
Xk.m- The sets C, and Z[i] denote the set of all complex 
numbers, and the Gaussian integer ring, respectively, where 
i = If z G C, then 5R(z) is the real part and is the 



imaginary part of z. Let | • | denote the modulus of a complex 
number. The || • || operation denotes the norm square of a 
vector. For a complex number z, the closest Gaussian integer 
to z is denoted by [z~\ , which we refer as the quantization of z. 
The notation [v] is the component-wise quantized version of 
the vector v. The Hermitian product of a and b is defined by 
(a, b) = h h a. Finally, the set of orthogonal vectors generated 
by the Gram-Schmidt orthogonalization procedure are denoted 
by {GS(d!),...,GS(d„)}. 

II. Background on Lattices and CLLL Algorithm 

A lattice A with basis {gi, g2, • • • , gn}, where gj, £ C™, is 
the set of all points {x = uG|u £ Z[z]™}. A generator matrix 
for A is an n x n complex matrix G = [gf, • ■ • , gn T ] • 
The Gram matrix of A is M = GG . The m-th successive 
minima of A, denoted by A m , is the radius of the smallest 
possible closed ball around origin containing m or more 
linearly independent lattice points. 

In complex lattice reduction, we let G' = UG, where U 
is an unimodular matrix. Let us define 

(g£,os ( g ;.)) 
m ~ IIGs( g ;.)ll 2 

where 1 < l,j < n. A generator matrix G' is said to be 
CLLL-reduced if the following two conditions are satisfied j2]: 

1) for 1 < j < £ < n 

< 1/2, < 1/2, 

2) for 1 < m < n, 

||GS (g' m )\\ 2 >(S-\^ n , m ^\ 2 ) ||GS (g^OH 2 

where 5 £ (1/2, 1] is a factor selected to achieve a good 

quality-complexity tradeoff. 
An algorithm is provided in J2) to evaluate a CLLL-reduced 
basis matrix G' of a lattice A with a generator matrix G. The 
input of this algorithm is the matrix G and a factor S, and the 
outputs of the algorithm are the unimodular matrix U and the 
CLLL-reduced basis matrix G' such that G' = UG. 

III. Signal Model and Problem Statement 

A flat-fading MIMO channel with n transmit antennas 
and n receive antennas is considered. The channel matrix H 
is in C nx ", where the entries of H are i.i.d. as CAf(0, 1) 
this channel coefficient remains fixed for a given interval (of 
at least N complex channel uses) and take an independent 
realization in the next interval. We use a n-layer transmission 
scheme where the information transmitted across different 
antennas are statistically independent. For 1 < m < n, the 
m-th layer is equipped with an encoder £ m : lZ k —> C N . This 
encoder maps a vector message m £ lZ k , where 7Z is a ring, 
into a lattice codeword x m £ A C C N . If X denotes the 
matrix of transmitted vectors, the received signal Y is given 
by Y = v^PHX + Z, where P = ^ and SNR denotes 
the average signal-to-noise ratio at each receive antenna. The 
entries of Z are i.i.d. and distributed as CAf(0, 1). We also 



assume that the channel state information is only available at 
the receiver. The goal of IF linear receiver is to approximate 
H with a non-singular integer matrix A. Since we suppose 
the information symbols to be in the ring 7Z, we look for an 
invertible matrix A over the ring 7Z. Thus, we have 

Y' = BY = VPAX + VP(BH A)X + BZ. 

As \/PAX is the useful signal component, the effective noise 
is \/ r P(BH — A)X + BZ. In particular, the power of the m-th 
row of effective noise is j|b,„|| 2 + P||b„,H — a m || 2 . Hence, 
we define 

5 (a m ,b m ) 4 ||b m i| 2 + P||b m H-a m || 2 , (1) 

where a m and b m denote the m-th row of A and B, 
respectively. Note that in order to increase the effective signal 
to noise ratio for each layer, the term g(a m ,h m ) has to be 
minimized for each m by appropriately selecting the matrices 
A and B. We formally put forth the problem statement below 

Given H and P, the problem is to find the matrices B £ 

C" x " and A £ Z\i] nxn such that 

• the maxi< m <„ <?(a m , b m ) is minimized, and 
« the matrix A is invertible over the ring 7Z. 

In order to have ZF receiver, we put B = H 1 and A = I„. 
If we let B = H^S" 1 where 

S = P~ l l n + HH ft , (2) 

and A = I„, then we get the linear MMSE receiver. 

In [141 . the authors have proposed a method to obtain A 
and B, which reduces g(a. m , b m ) for each m. We now recall 
the approach presented in lfl4l . First, conditioned on a fixed 
a m = a, the term g(a, b m ) is minimized over all possible 
values of b m . As a result, the optimum value of b m can be 
obtained as 

b m = aH^S" 1 . (3) 

Then, after replacing h rn of <£3j in g(a, b rra ), one has to 
minimize g(a, aH^S -1 ) over all possible values of a to obtain 
a m as a m = argmin a g(a, aH^S -1 ). The last expression can 
be written as 

a m = argmin aVDV^a' 1 , (4) 

a 

where V is the matrix composed of the eigenvectors of 
HH' 1 and D is a diagonal matrix with m-th entry D m m = 

(PPm + l) ' where Pm is the m-th singular value of H. 
With this, we have to obtain n vectors {a m } which result in 
the first n smaller values of aVDV^a' 1 along with the non- 
singular property on A. In order to get a m , 1 < m < n, the 
authors of lfl4l have suggested an exhaustive search for each 
component of a m within a sphere of squared radius 

1 + Ppl a » (5) 

where p max = max m p m . It has also been pointed out in lfT4l 
that this search can be accelerated by means of a sphere 
decoder on the lattice with Gram matrix M' = VDV' 1 . For a 



fixed P, the complexity of this approach is of order O (P n ). 
It is also shown in [14] that the exhaustive search approach 
provides a diversity order of n and a multiplexing gain of n. 
The outage rate and probability of this scheme are studied 
in fl4l . At this stage, we note that the exhaustive computation 
of a m has high computational complexity, especially for large 
values of P and n, and hence the approach in lfT4l is not 
practical even for the 2x2 complex case. 

The best possible a m 's for this problem are the set of all 
successive minimas of VD2, which can be approximately 
computed by either the Minkowski or HKZ lattice reduction 
algorithms Q, 0, 02). In 0, using HKZ and Minkowski 
reduction algorithms, we obtain the matrices A and B and 
show that these practical algorithms achieve full receive diver- 
sity in terms of error performance with reasonable complexity. 
Since the CLLL algorithm has lower complexity than HKZ 
and Minkowski reduction algorithms, it could be employed to 
get matrix A as well. Hence, in the rest of the paper we focus 
on using CLLL algorithm to find the best possible A for our 
problem. 

IV. Low-complexity IF Receivers 

In this section, we propose three low-complexity methods 
to obtain some candidates for the rows of A. The first two 
are based on CLLL algorithm and the last one is based on 
SVD decomposition. Then, we propose a selective combining 
technique to choose the rows of A from the candidate rows. 
Once we obtain A, we obtain B as B = AH^S -1 , where S is 
given in (0. Henceforth, we only address systematic methods 
to obtain A. 

A. Algorithm 1 via CLLL 

It is pointed out in lfl4l that the minimization problem in 
is the shortest vector problem for a lattice with Gram matrix 
M' = VDV' 1 . Since M' is a positive definite matrix, we can 
write M' = LL h for some L G C«x" by using Choelsky 
decomposition. With this, the rows of L = VD2 generate 
a lattice, say A'. Based on (@), a set of possible choices 
for {ai, . . . , a„} is the set of complex integer vectors whose 
corresponding lattice points in A' have lengths at most equal 
to the n-th successive minima of A'. However, finding these 

Input: H G C" x ™, and P. 

Output: A set of n candidates for the rows of A. 

1) Form L = VD^ of a lattice A'. 

2) Reduce L to L cUl by CLLL algorithm. 

3) The n rows of L ciH L _1 provides n values of a m . 

We now use the CLLL algorithm |2] to obtain the complex 
integer vectors. In particular, we use the CLLL algorithm to 
reduce the basis set in L to obtain a new basis set represented 
by the rows of h clu . For each 1 < m < n, it is known that 
the length of the m-th row vector in L ciH is upper bounded 
by a scaled version of the m-th successive minima of A' 0. 
Hence, the rows of L ciH L _1 can be used to obtain n possible 
choices for the desired matrix A. We note that the structure of 
the above algorithm is exactly the same as the one presented 



in 00 with only one difference: we used HKZ and Minkowski 
lattice reduction algorithms instead of CLLL algorithm. 

B. Algorithm 2 via CLLL 

Given H and P, let us define a 2n— dimensional complex 
lattice A generated by 
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The Gram matrix of this lattice is M = GG''. The Schur's 
component of M denoted by (M|S) is given by (M|S) = 
I n — H' l S" 1 H, while the Schur's decomposition of M can be 
derived as, 
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Replacing H by its SVD representation, one can easily observe 
that M' = (M|S). Consider x G A where x = uG and 
u G 1[i] 2n and u = [b|a] forming by adjoining row vector 
a after b for some a,b G Z[i] n . For u = [d|c] with 
d, c G Z[i] n , we define /(c, d) = uMu'\ Further, we expand 
the term /(c,d) to obtain /(c,d) = P -1 ||d|| 2 + ||dH - c|| 2 . 
Note that /(c, d) = P~ l g(c, d). Therefore, the solutions 
to the minimization of f(c, d) are also the solutions to the 
minimization of g(c, d). Towards minimizing f(c, d), one can 
immediately recognize that solving minimization of /(c, d) is 
nothing but finding the shortest vector of A. For the matrix 
G in ©, let Q clu denote the 2n-dimensional CLLL -reduced 
generator matrix of A. Using the short vectors in G cUl , we 
obtain the complex integer matrix U with 2n row vectors 
u m = [dm I Cm] such that U = G cHi G _1 . With this, we 
have obtained complex integer vectors d m , c TO which results 
in smaller values for g(d, c). Hence, the vectors d m ,c m can 
be readily used for the IF architecture as h m = d m , and 
a m = c m - However, we do not use h m as it is a complex 
integer vector. In such a case, b m H ^ a m even for large SNR 
values which in turn results in an error floor in the probability 
of error. Instead, we only use a,„, and subsequently obtain b m 

Input: H G C" x ", and P. 

Output: A set of 2n candidates for the rows of A. 

1) Form the matrix G as in (0. 

2) Reduce G to G c/ " by CLLL and compute U = 

3) Use u m = [d m I c OT ], choose a m = c m for 1 < 
m < 2n and put b m as in ([3]). 

Remark 1: It is important to highlight the difference be- 
tween Algorithm 1 and Algorithm 2. In Algorithm 1, a set 
of candidates for a m are chosen via the CLLL reduction 
technique on a n-dimensional complex lattice. However, in Al- 
gorithm 2, we jointly obtain d m and c TO on a 2n-dimensional 
complex lattice, and only use a m = c TO as candidates for the 
rows of A. Note that vectors a m delivered from Algorithm 2 
can be different from that of Algorithm 1, since Algorithm 2 
attempts to minimize g(a, b), whereas Algorithm 1 attempts 
to solve minimization of aVDV' l a'\ 



Note that, the CLLL algorithm guarantees the existence of at 
least n vectors of the lattice A such that the vector in the first 
row of R is most likely the shortest vector of the lattice. 

C. Algorithm 3 via SVD 

For large values of P, we can |[T4| write a m VDV' l aJ^ in 

© as p\ 2 ||via^|| H h p~ 2 ||v„a^|| , where v fc denotes 

the fc-th row of V. From the SVD property, we have p m i n = 
pi and p m ax = p n - If Pmax > Pj for j ^ 7i, then the above 
equation suggests us to select all a m , 1 < m < n, along the 
direction of v„, and as short as possible. On the other hand, 
if pk is large for each k, and are comparable, then a set of 
good candidates can come from choosing a complex integer 

Input: H e C" x ™, and P. 

Output: A set of n candidates for the rows of A. 

1) Obtain the SVD of H as H = VEV h . 

2) Choose a m = |_ v m] f° r 1 < m < n. 

Therefore, another n possible choices for integer vectors can 
be obtained from the rows of V as a TO = [ v ml f° r 1 < m < 
n, where [•] denotes the nearest integer of a real number. 

D. Combined CLLL-SVD solution 

Till now, we have proposed three different low-complexity 
algorithms to obtain candidate vectors for the rows of A. 
Algorithm 1 gives us n choices, Algorithm 2 delivers 2n 
choices, whereas Algorithm 3 brings another n possible 
choices. Overall, we can use all the An candidate vectors for 
a m and obtain the corresponding vectors b r „, 1 < m < An 
using the relation in ([3J. With this, we have An pairs of 
(a m ,b m ). Now, we proceed to sort these An vectors in the 
increasing order of their g(a m , b m ) values. Then, we find the 
first n vectors which form an invertible matrix over the ring 
1Z. Note that all the three algorithms have low computational 
complexity. As a result, the combining algorithm has lower 
computational complexity than the exhaustive search and other 
lattice reduction algorithms like HKZ and Minkowski, and 
hence, the proposed technique is amenable to implementation. 
We refer to the combined solution as "combined CLLL-SVD 
solution". 

E. Complexity Comparison 

It is clear that the combined CLLL-SVD method includes at 
most a SVD algorithm and a matrix inversion with complexity 
0(n 3 ), a 2n-dimensional CLLL algorithm with complexity 
O [An 2 log(2n)) following by a sorting algorithm of size An 
with complexity O (4nlog(4n)). Therefore the complexity of 
this approach is O (n 3 ) and it is obviously independent of 
P. The complexity of the proposed exhaustive search in [fT31 
is known to be 0(P n ). The complexity of other reduction 
algorithms like HKZ and Minkowski are given in (8j. In 
Table H we compare the complexity of finding matrix A for 
different known methods. 
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Fig. 1 . Ergodic rate of various linear receivers for 2x2 MIMO channel. 



V. Simulation Results 

In this section, we present simulation results on the ergodic 
rate (see (8) or J9) for definition of ergodic rate) and the 
probability of error for the following receiver architectures on 
2x2 MIMO channel: (/) IF linear receiver with exhaustive 
search, (ii) IF linear receiver with Minkowski lattice reduction 
solution (8), (Hi) IF linear receiver with combined CLLL-SVD 
solution, 0'v) the ZF and MMSE linear receiver, and (v) the 
joint maximum likelihood (ML) decoder. For the IF receiver 
with exhaustive search, the results are presented with the 
constraint of fixed radius for the exhaustive search. We have 
not used the radius constraint given in (0 as the corresponding 
search space increases with P. Instead, we have used a fixed 
radius of 8 for all values of P. By relaxing this constraint, we 
have reduced the complexity of brute force search, noticeably. 

In Fig. [U we present the ergodic rate of the above listed 
receivers, wherein, for the case of ML receiver, the ergodic 
capacity |[T2l of 2 x 2 MIMO channel has been presented. The 
rate for IF linear receiver with Minkowski lattice reduction 
solution is cited from J5). We observe that the combined 
CLLL-SVD solution performs pretty much the same as IF 
receiver with exhaustive search and IF based on Minkowski 
lattice reduction solution at low and moderate SNRs. Also, 
note that the combined CLLL-SVD and Minkowski lattice 
reduction solutions give lower bounds on the ergodic rate of 
exhaustive search based IF receiver while the latter one is 
tighter. For the ergodic rate results of the IF receiver, we have 
used A matrices which are invertible over Z[i]. 

Now, we present the uncoded bit error rate (BER) for the 
above receiver architectures with 4-QAM constellation. We 
use the finite constellation S = {0, 1, i, 1 + i} carved out of 
the infinite lattice Z[i]. Note that Z[i] = S © 2Z[z], In this 

/ SNR 

method, the received vector is of the form, y = \J ^Hs+z, 
where s £ S 2xl . For the above setting, we use modulo lattice 
decoding at the receiver. Each component of y is decoded 
to the nearest point in Z[i] and then "modulo 2" operation 
is performed independently on its in-phase and quadrature 
component. With this, we get r £ S 2 x 1 from both the com- 



TABLE I 
Complexity Comparison. 



Name 


Approach 


Complexity 


Integer-Forcing Linear receiver 


Brute Force 
HKZ reduction 
Minkowski reduction 
Combined CLLL-SVD solution 


pn 

(27re) ,l+ °( log2 ") 
(5/4) 2 ™ 2 
n 3 




ML 

IF (fixed radius) 

IF Minkowski 

IF Combined CLLL-SVD 

MMSE 

ZF 
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Fig. 2. BER for various lineal' receivers with 4-QAM constellation. 



ponents of y. Further, we solve the system of linear equations 
r = As modulo 2 over the ring 7Z to obtain the decoded vector 
s. In Fig. [2] we present the BER results for all the five receiver 
architectures. Note that both the IF receiver with combined 
CLLL-SVD and Minkowski lattice reduction solutions out- 
perform the ZF and MMSE architectures, but trades-off error 
performance for complexity in comparison with brute force 
search. In particular, the combined CLLL-SVD approach fails 
to provide diversity results as that of the exhaustive search and 
Minkowski lattice reduction approaches. This diversity loss 
is due to the larger value of max m g(a m , b m ) delivered by 
the combined CLLL-SVD algorithm in comparison with the 
optimum solution. For the probability of error results of the 
IF receiver, we have used A matrices which are invertible 
over %2 [i] ■ We have observed similar results for 4x4 MIMO 
channels as well. 

VI. Conclusions 

In |8j, Algorithm 1 along with HKZ and Minkowski lattice 
reduction algorithms was employed to find the matrix A. 
Algorithm 1 which includes the complex LLL method as an 
alternate technique for the above reduction methods turned 
out to be not satisfactory in terms of ergodic rate and error 
performance. Hence, we have tried to improve effectiveness 
of CLLL algorithm by combining it with other approaches 
such as Algorithms 2 and 3. We have proposed a low- 
complexity systematic method called the combined CLLL- 
SVD algorithm for the MIMO IF architecture. Simulation 



results on the ergodic rate and the probability of error were 
also presented to reveal the effectiveness of combined CLLL- 
SVD solution versus other linear receivers. The proposed com- 
bined algorithm trades-off error performance for complexity 
in comparison with both IF receivers based on exhaustive 
search and Minkowski or HKZ lattice reduction algorithms. 
Further improvements are required to achieve results which 
are competitive with IF receivers based on exhaustive search 
and the ones presented in 0. 
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